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Abstract 



^ ' LABSMC Monte Carlo event generator is used to simulate Bhabha scattering at 

^D . high energies. Different sources of radiative corrections are considered. The resulting 

ji^ I precision is discussed. 



PACS: 12.20.-m Quantum electrodynamics, 12.20.Ds Specific calculations 



o 

_ 1 Introduction 

p ' In this note we are going to discuss the application the Monte Carlo event generator LABSMC [Q 
to the case of LEP2 energies. The topic is actual now in view of the analysis of LEP2 data. 
The experiment requires high precision theoretical predictions to perform more deep tests of 
the Standard Model and to look for a new physics. 
r> I Initially the LABSMC event generator was developed to simulate large-angle Bhabha scatter- 

ed ■ ing at energies of about a few GeV's at electron positron colliders like VEPP-2M and DA$NE. 
The code included the Born level matrix element, the complete set of 0{a) QED RC, and the 
higher order leading logarithmic RC by means of the electron structure functions. The rele- 
vant set of formulae can be found in Ref. 0. The generation of events is performed using an 
original algorithm, which combines advantages of semi-analytical programs and Monte Carlo 
generators. 

The structure of our event generator was described recently in paper [|l|. The extension for 
higher energies is done by introducing electroweak (EW) contributions, such as Z-exchange, 
into the matrix elements. The third [0] and fourth |^ order leading logarithmic photonic 
corrections were also included in the new version. So, the structure of the code is kept the 
same, and we have to describe now what kind of EW effects were included in our code to 
work with large-angle Bhabha at LEP2. In particular we are going to consider the region 
of radiative return to the Z-peak in radiative Bhabha scattering. To estimate the resulting 
theoretical uncertainty, one has to analyse various sources of radiative corrections (RC). 

*on leave of absence from Joint Institute for Nuclear Research, Dubna, Russia. 



2 Electroweak contributions 

The set of electroweak effects was included according to Refs. P, P, 0. The Born level cross 
section in the program contains now both the photon and Z-boson exchange contributions (it 
is used also as a kernel cross section for higher order leading log radiative corrections). The first 
order virtual and soft EW RC were taken directly from the semi-analytical code ALIBABA |^. 
The EW matrix element for radiative Bhabha is taken from Ref. 0. A comparison of the 
Bhabha cross-section integrated over photons is in a reasonably good agreement (see Table 1, 
error-bars are dropped) with the published numbers of other codes p]. 

Table 1: Comparison with Fig. 21 from Ref. p[, cross-sections in pb. 



EcM, GeV 


BHWIDE 


TOPAZO 


BHAGENE3 


UNIBAB 


SABSPV 


BHAGEN95 


LABSMC 




^acol = 10° 


175 
190 
205 


35.257 
29.899 
25.593 


35.455 
30.024 
25.738 


34.690 

28.780 
24.690 


34.498 
29.189 
25.976 


35.740 
30.270 
25.960 


35.847 
30.352 
26.007 


35.337 
29.941 
25.687 




^acol = 25° 


175 
190 
205 


39.741 
33.698 
28.929 


40.487 
34.336 
29.460 


39.170 
32.400 

27.840 


39.521 
33.512 

28.710 


40.240 
34.100 
29.280 


40.505 
34.331 
29.437 


40.029 
33.954 

29.178 



3 Radiative return with a visible photon 

At LEP2 the radiative return to the Z-peak due to photon or pair radiation gives a sizable 
contribution to the cross-section. This process is used itself in particular to look for anomalous 
gauge boson couplings. 

The pure tree level matrix element |^ for radiative process 
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+ e' 
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was supplemented by radiative corrections due to initial state soft and hard collinear radiation 
by means of the electron structure function approach . The electron-positron pair produc- 
tion was taken into account in the same way. As an energy scale for the structure functions we 
took the t-channel momentum transferred, because the corresponding diagrams are dominant. 
The vacuum polarization correction to the photon propagators is apphed as well. 

In Table 2 we put the result for the following conditions: ii^cM=183, 189 GeV; | cos 6*6*1 < 
0.95; at least one electron has | cos^el < 0.7; electrons should have transverse momenta above 
1 GeV; the final particles are to be isolated by at least 20 degree from each other; the total 
observed energy > 0.8 Eqm] | cos 6*^1 < 0.7. In the columns without Z-peak we excluded 
the events with the invariant mass of the electron positron pair in the range 85 GeV < 
Mee < 95 GeV. As could be seen from the numbers, the ISR LLA corrections are in this case 
of the order 2%. The additional non-standard LLA corrections, which were found in Ref. [B!Ol 



make a small shift of the correction; but an independent verification of the investigation is 
required. 



The complete set of 0{a) EW radiative corrections to the process (|1|) is unknown. To 
estimate the uncertainty of our result we look at the relative size of the known leading and 
sub-leading 0{a) RC to the Bhabha process itself. For an analogous set of cuts for Bhabha 
scattering for the difference of the correction values we have 6tot — (^lla ~ 1%Q ^^ this way 
we estimate the precision of our results for the radiative process (|ip to be of the order 1.5%. 



Table 2: The cross section in pb of radiative Bhabha with a visible photon in different ap- 
proximations. 





total 


without Z-peak 


EcM GeV] 


183 


189 


183 


189 


tree-level 

vacuum polarisation 

vac. pol. + ISR LLA 


0.9817 
1.1022 
1.0842 


0.9146 
1.0342 
1.0088 


0.8251 
0.9630 
0.9346 


0.7727 
0.8853 
0.8770 



4 Conclusions 



The precision of the theoretical predictions, which can be obtained by means of the presented 
code for inclusive large-angle Bhabha scattering at LEP2, is estimated to be of the order 0.2%. 
It is defined mainly by the unknown radiative corrections of the order 0{{a/7r)'^L) ?» 10~^. 
The coefficients before these terms are not too large, as was seen in the case of small-angle 
Bhabha scattering. At LEP2 energies in the large-angle Bhabha process the t-channel photon 
exchange is dominant. Nevertheless, a correct treatment of the electroweak Born and the first 
order corrections is important. The technical precision of the code is to be verified in further 
test and comparisons with other codes. By means of the comparisons of the semi-analytical 
branch and the pure Monte Carlo one of the code we have a good control of such parameters 
as the precision of numerical integration and the number of events to be generated. That 
allows to reach an ordered level of the uncertainty in numerical evaluations. 

Another source of uncertainties is an incomplete treatment of pair production. In the 
current version of the code the pair production corrections are evaluated in the 0{a'^L'^) 
structure function approach. Both the singlet and non-singlet electron pairs are included. A 
special brunch of the code to scrutinise the pair production |]TT| is in progress. 

The inclusion of the third and fourth order LLA photonic corrections allows not to use 
exponentiation. A simple estimate shows that the difference between the two treatments 
at LEP2 is negligible, while the exponentiation requires a quite different event generation 
procedure. 
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^The special cut on the scattering angle of "at least one electron has |cos6'e| < 0.7" is similar to the 
narrow-wide event selection in small-angle Bhabha at LEPl. In both cases we see a considerable reduction 
of the RC size. If we apply this cut, the difference ^tot — (^lla is well below the 1% level. 
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